Occupation time fluctuation limits of infinite 
variance equilibrium branching systems 



We establish limit theorems for the fluctuations of the rescaled occupa- 
tion time of a (d, a, /3)-branching particle system. It consists of particles 
moving according to a symmetric a-stable motion in R d . The branch- 
ing law is in the domain of attraction of a (l+/3)-stable law and the 
initial condition is an equilibrium random measure for the system (de- 
fined below). In the paper we treat separately the cases of intermediate 
a/f) < d < (1 + f3)a/f3, critical d = (1 + P)a/P and large d > (1 + (3)a/P 
dimensions. In the most interesting case of intermediate dimensions we 
obtain a version of a fractional stable motion. The long-range depen- 
dence structure of this process is also studied. Contrary to this case, limit 
processes in critical and large dimensions have independent increments. 
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1 Introduction 

1.1 Branching system and occupation time fluctuations 

The aim of this paper is to present some (functional) limit theorems for the 
occupation time fluctuation process of a branching particle system. We call a 
(d, a, /^-branching particle system (denoted in the sequel by N) a set of particles 
moving independently according to the spherically symmetric a-stable Levy 
motion (0 < a < 2) in R d and splitting after exponential time (with intensity 
V) with branching law 
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Abstract 




(0 < (3 < 1). This is an example of a law in the domain of attraction of (1 + /in- 
stable variable. It has infinite variance and is critical. For f3 — 1 it reduces 



to binary critical branching which was treated in a series of papers mentioned 
below. The generating function of this law is 



1 + p 

The particle system will be represented by an empirical measure process {N t )t>o, 
i.e. for a Borel set A, N t (A) is a (random) number of particles in A at time t. 
The initial particle distribution is yet to be introduced. The most natural choice 
is a Poisson random field with homogeneous intensity, i.e. (Lebesgue measure) 
A. This case, which was studied in 0) and (d), is a starting point and reference 
for our investigation. It is known (|13T ) that for a/(3 < d such system (denoted 
by 7V Polss ) converges to an equilibrium distribution 

N Pmss Eq (1.2) 

where => denotes weak convergence in the space of point measures. The Laplace 
functional of the equilibrium distribution is given by 



Eexp {- {Eq, <p)} = exp U \, e-«> - 1) + V J o (X,H (j (-,*))) ds\ , (1.3) 
where 

j(x,l):=Eexp(-{N?,<p)) (1.4) 

N x is the empirical process of the system starting from x G R d , H(s) — F(s) — s, 
if : R d -> R+, ip e C 1 ^) n C(R d ) and j satisfies the integral equation 

j (x, I) = (x) + V [ T t _ s H (j (•, s)) (x) ds. 



This equations can be obtained in the same way as (0, (2.4)). In this paper we 
consider a system N starting off from Eq and compare the obtained result to the 
ones in (0) and H). For the process (N t )t>o we define the rescaled occupation 
time fluctuations process by 

1 f Tt 

X T {t) = — (N s -X)ds, (1.5) 

b T Jo 

where Ft is a proper normalization and T is a scaling parameter which accel- 
erates the time. The object of our investigation is the limit of Xt as T tends 
to +oo 

X T => X. (1.6) 

For the time being we are not very rigorous and do not specify the type of 
convergence. 
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1.2 Results and proof techniques 

In the proofs we will rely on methods presented in H), (0) and (flih. 
Although the process X? is signed-measure-valued it is convenient to regard it 
as a process with values in the space 6>'(IR d ) of tempered distributions, which is 
dual to the space of smooth and rapidly decreasing functions <S(R d ). We denote 
duality in this space by (•, •}. In this space one may employ space-time method 
introduced by which together with Mitoma's theorem constitute a powerful 
technique in proving weak, functional convergence. 

Three kinds of convergence are used. The convergence of finite-dimensional 
distributions is denoted by For a continuous, S"(R d )-valued process X — 
(X t )t>o and any r > one can define an 5'(R d+1 )-valued random variable 

(x,*) = jf {X s ,$(-,s))ds, $e5(M d+1 ). (1.7) 

If for any r > X n — > X in distribution, we say that the convergence in 
the space-time sense holds and denote this fact by =>-j. Finally, we consider 
the functional weak convergence denoted by X n => c X. It holds if for any 
r > processes X n — (X n (t))t e [o jT ] converge to X = (X (t)) te [o tT ] weakly in 
C([0, t], 5'(R d )) (in the sequel without loss of generality we assume r = 1). It 
is known that =>j and =>/ do not imply each other, but either of them together 
with tightness implies => c Qj. Conversely, =4> c implies both =>j, =^f. 
The presentation of the results naturally splits into parts, corresponding to 
intermediate dimensions 



critical 



and large dimensions 



a , a(l+/3) 

' 
d> a(l+0) 



respectively. 

In the first case of intermediate dimensions we obtain a weak functional con- 
vergence to process X of the form X = Kr/X where if is a constant, r/-a stable 
process being in a sense a stable (non- Gaussian) analogue of a fractional Brow- 
nian motion. So we see that in this case the limit has a very simple spatial 
structure whereas its temporal structure is complicated. It is also worthwhile to 
point out that 77 has a stationary increments (unlike the corresponding process 
in (0)) and is a heavy-tailed process with long-range dependence. This depen- 
dence is described in Section [2] in terms of dependence exponent and roughly 
speaking it means that " dependence" decays polynomially. The cases of critical 
and large dimensions differ substantially, one can prove only finite-dimensional 
distributions and space-time convergences. In both cases we obtain processes 
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with independent increments and the limit for large dimensions is truly S' (Re- 
valued. These processes are not continuous hence it is not possible to obtain 
functional convergence. 

1.3 A survey of results 

This paper is a part of a larger programme carried out by Bojdecki et al. and 
recently by Milos. It seems useful to present it here in a compact way. This 
small survey is not meant to be exhaustive nor very strict it aims only to give 
a reader a glimpse of the whole picture. In tables below we gather limits of 
occupation time fluctuations under time rescaling (as defined by (|1.6[) ). type 
of convergence and the normalizing factor in different settings. The structure 
of tables reflects dependence on the dimension of the space and the starting 
distribution. 



Table 1. Systems with finite variance branching law 





Poisson 


Equilibrium 


a < d < 2a 
intermediate 


K ■ sub-frac-BM- A 

m 3-d/a 

functional T 2 
(5) and (16) 


K ■ frac-BM ■ A 

m 3 - d / a 

functional T 2 
(16) 


d = 2a 
critical 


K ■ BM ■ A 
functional (T log T) i 
(6) and (15) 


the same 
(15) 


d>2a 
large 


S'{R d )-BM 
functional T2 
(6J and QJ5) 


the same* 
05) 



* - due to technical difficulties functional convergence proved only for finite 

fourth-moment branching law. 

K denotes a generic constant 

BM - standard real-valued Brownian motion 

sub-frac-BM - sub-fractional Brownian motion (ie. centered Gaussian process 
with covariance function s h + t h - |[(s + t) h + \a — t\ h ]) 

frac-BM - fractional Brownian motion (ie. centered Gaussian process with co- 
variance function ^[s h + t h + \s — t\ h ]) 

S' (M. d )-BM - centered Gaussian 5'(IR d )-valued process with covariance func- 
tional 



Cov ((X s , cpi) , (X it w)) = (s A *)— + ^pj <pi{z) V2 {z)dz, 

where <fi,if2 G 5 (K d ) and m depends on branching law. Papers |j§) and (@) 
contain also results for systems without branching. 
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Table 2. Systems with infinite variance branching law - generating function (jl.lj) 





Poisson 


Equilibrium 




K ■ sub-frac-SM- X 


if • frac-SM ■ X 


a „ j „ a(l+/3) 

intermediate 


functional 


functional 


F T = r(2+/3-#/3)/(i+/3) 






(7) 


this paper 




K-SM- X 




a— p 
critical 


fin-dims and space-time 


the same 


(TlogT)TTF 

(S) 


this paper 




5'(M <t )-5M 




a > ^ 
large 


fin-dims and space-time 


the same 


(8) 


this paper 



Here, 

sub-frac-SM - "sub-fractional" stable motion, defined by (|2.2[) 
frac-SM - "fractional" stable motion, defined by ()2.4p 

SM - stable motion with independent increments, with finite dimensional dis- 
tributions given by (|2.5j) 

S' (M. d )-SM - 5'(IR !i )-valued stable motion with finite dimensional distributions 
given by (|2.6p . 

Let us notice first that the results for the case of finite and infinite variance are 
in a sense similar. The processes in Table 1 are Gaussian counterparts of stable 
processes in Table 2. Informally speaking, the finite variance branching law is 
a limit of laws given by (jl.ip hence one can observe similar phenomena in both 
cases. The case of intermediate dimensions is most interesting. The limits have 
similar spatial structure and complicated temporal one with long-range depen- 
dence property. The dependence on the starting distribution is intriguing since 
Eq measure is the limit for a Poisson-starting system (jl.2p and by the time did 
not acquire any intuitive explanation. However both limits have the long-range 
dependence property but for the equilibrium-starting system this dependence is 
stronger. 

The remarkable feature of the limit process in the equilibrium case is that it can 
be decomposed into a sum of two independent (1 + /9)-stable processes. One of 
them being exactly the limit process in the Poisson case. This decomposition is 
an analogue of the one studied in (|12| ) for fractional Brownian motion. It should 
be noted that the process obtained in this case is a stable analogue of fractional 
Brownian motion. Namely, it is self-similar and has stationary increments. Pro- 
cesses with this properties were discussed in (|18l . Chapter 7) (see also Remark 
12.71) . This fact makes analogies between infinite- variance and finite- variance 
cases even stronger (recall that the limit process in the case of the finite vari- 
ance branching law is fractional Brownian motion). The cases of critical and 
large dimensions are less complicated. The limits have independent increments 
and complicated spatial structure for large dimensions going beyond the space 
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of measures. The qualitative change of the type of limits with dimension can 
be, partially, explained by recurrence and transient property of the underlying 
a-stable Levy motion. Further results on the fluctuations of the occupation 
time can be found in (11), (fioh . (0) were high-density limits and system with 
inhomogeneous starting distributions are studied. One should also mention (3) 
and (0) where similar problems are considered in discrete setting (lattice Z d ). 
Although the proofs in this paper rely mostly on the schema and methods used 
m (0), (H) and (0) we had to overcome some new technical difficulties which 
emerged during studies on additional terms arising in analysis of equilibrium- 
starting system. 



2 Results 

By T we denote the semigroup of the a-stable motion and by pt its transition 
density, i.e., 

T t f(x) = (p t *f)(x). (2.1) 

Let M be an independently scattered random (1 + /3)-stable measure on 1^x1+ 
with the Lebesgue control measure. More precisely, for a Borel set A, M(A) is 
(1 + /3)-stable variable with characteristic function 

explA^lzl 1 ^ (l - isgn(z)tan|(l +£))} , 

variables on disjoint set are independent and M is er-additive a.s. 
We define two stable processes. 

Vt= [ (l[o,t](r) [ Pu-r(x)du] M(dx,dr) (2.2) 
jR d +i V J r J 

and + ^ 

^= I 1(1 Ps+i(x)ds) M(dx,dl) (2.3) 

JO JR d \Jo / 

It can be checked that for intermediate dimensions both processes are well- 
defined in the sense given in d Chapter 3). 
Assume now that rj 1 and rf are independent, then we define 

T) := r] 1 + rf . (2.4) 

This process plays fundamental role in this paper. Detailed presentation of its 
properties is postponed to Theorem l2.6l and Proposition ^. 91 Now we give series 
of three theorems which are the main results of the paper. 
In these theorems Xt is the rescaled occupation time fluctuation process defined 
by (|1.5|) for a system N starting from equilibrium distribution (|1.3[) . 

Theorem 2.1. Assume | < d < and F T = r( 2 +' 3 -^)/( 1 +' 3 ) . Then 
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wher 



V 7T X 1+/3 



Remark 2.2. As announced in the Introduction the process n consists of two 
independent summands (|2.4J) where n 1 is the process that occurs in the limit 
for the Poisson system (see ([7j)). An explanation of the reason of this structure 
of rj as well as the interpretation of rf require further studies. If (3 = 1 and 77 
is a fractional Brownian motion on the the whole line then r\ x = ??t+ 2 ''~ t (sub- 
fractional Brownian motion) and 7] 1 = ' 7t ~ 2 ' ? ~ t (see (fl2h ). 

Theorem 2.3. Assume d = and F T = (TlogT)TTF. Then 

Xt —>i K\£ and Xt — ►/ -ftTA£ as T — > +00 

where £ is (1 + (3) -stable process with stationary independent increments and 
characteristic function 

Eexp(iz£ t ) = exp j-t|z| 1+/3 (l - isgn(z) tan |(1 + /?)) | , z G R, f > (2.5) 

and 

K=(-Vcoa^(l + 0)J^J o p r (x)dr^j Pl {x)dx^j 

Before presenting the last theorem we introduce the potential operator cor- 
responding to the a-stable motion. 



r+oo 

g/(x) = / T t f{x)dt. 

Jo 



Theorem 2.4. Assume d > a(1 + ffl and F T = TW. Then 

Xt — H X and Xt — >/ X as T — > +00, 

where X is an S' (M. d ) -valued (1 + f3)-stable process with stationary independent 
increments and characteristic function 

Eexp(i (X(t), (/))) = exp j-X 1+/3 i J \G<f)(x)\ 1+l3 (l - i{sgng<t>{x) tan |(1 + /?))) dxj , 

G <S(M d ),t > 0,(2.6) 

K ={-—r-2 {1+i,) ) 

Remark 2.5. The limits in the last two theorems have independent increments 
and are non-Gaussian hence by (|l4 , Theorem 13.4) are not continuous. It is 
somehow unexpected since processes Xt are clearly continuous. This is also the 
reason why we can not obtain functional convergence in those cases. 
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Propositions below summarize basic properties of 77 denned by (|2.4[) . 

Theorem 2.6. r\ is a continuous, (1 + @) -stable process process. It is self- 
similar with exponent H = (2 + [3 — ^)/(l + /3) and has stationary increments. 

The self-similarity can be proved by a simple calculation using the character- 
istic function of the finite-dimensional distributions of (12. 4p obtained using (l8l . 
(3.2.2)). Processes have stationary increments which comes straightforward 
from the fact that N is stationary (since it is a Markov process starting from 
stationary distribution) . From Theorem 12.11 we know that the process r\ is a 
limit of Xt hence is also stationary. 

Remark 2.7. Processes of this type were discussed in fl8l . Chapter 7). In the 
notation used there 77 is -ff-ssi stable process. Contrary to the Gaussian case, 
where there is a unique up to a constant ii-ssi process for a given H (fractional 
Brownian motion with Hurst parameter H), there are plenty of stable ii-ssi's. 
It would be interesting to check if 77 is one of already known processes (this 
could draw analogies to other problems) or is a new process. Unfortunately we 
do not know the answer to this question. 

We introduce now a general notation to investigate long-range dependence 
in the case of stable processes (see (0)). 

Definition 2.8. Let r\ be a real infinitely divisible process. For < u < v < 
s < t, T > 0, z%, Z2 eR define 

D T {z 1 ,Z 2 ;U,V,S,t) = \\ og ^i(^-ri^z 2 (r 1T+t -r 1TJts ) 

-logEe™ 1 ^-^) _ Ee «2(r,T +t -^T +s )| j (2.7) 
Dependence exponent n is defined by 

k = inf inf sup{7 > : Dt(zi, Z2, u, v, s, t) = o(T -7 ) as T — » +00} 
By (0, Theorem 2.7) exponent k of 771 (denoted by £ therein) is 

) P<lk- 

Conducting similar computations as in (7, Proof of Theorem 2.7) it can be 
checked that dependence exponent for 772 is k = ^ — 1. It is straightforward 
consequence of Definition 12.81 that the dependence exponent of a sum of in- 
dependent processes is minimum of the exponents of the summands hence we 
obtain 

Proposition 2.9. Process 77 has dependence exponent k = — — 1. 

Remark 2.10. It is interesting to notice that addition of an independent term 
772 arising in the limit for the equilibrium-starting system (recall Theorem 12. II) 
increases long-range dependence and the dependence exponent does not depend 
on (3 any more. 
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3 Proofs 



For the sake of brevity proofs for Theorem 12.31 and 12.41 are omitted. They are 
direct combination of the methods of (|8|) and the argument employed in the 
proof of Theorem 12.11 The scheme below is quite general and could be easily 
adapted for those proofs. 



3.1 Scheme of the proof 

To make the proof clearer we present a general scheme. Detailed calculation are 
deferred to a separated section. We treat measure- valued processes as £>' (Re- 
valued one. This enables usage of space-time method from Q). Let Xt denote 
5'(M d )-random variable defined by (|1.7[) corresponding to the process Xt- In 
order to prove weak convergence in C([0, l],5'(M d )) to X it suffices to prove 
weak convergence of Xt 

X t ^X (3.1) 
and tightness of {Xt}t>i- In order to obtain (|3.ip it suffices to verify that 

Ee <*r,*) ^ Ee (*,*) (3.2) 

for any non- negative <E> € S(R d ). The tightness can be proven utilizing the Mit- 
oma theorem (jilt ), which states that tightness of {Xt}t>i m C([0, T],<S'(R <i )) 
is equivalent to tightness of (X T , 4>) in C([0, r],R d ) for any 4> G S(R d ). 

3.1.1 Space-time convergence 

The purpose of this subsection is a calculation of the Laplace transform and 
gathering facts used to show convergence (|3.2p . The schema described below 
generally follows the lines of a scheme presented in (0), Q) and (flih hence we 
omit some details. 

To make the proof shorter we will consider $ of the special form: 

$(x, t) = <p(z)i>(t) <p £ S(R d ), tp e 5(R+), Lp > 0, > 0. 
We also denote 

Vt = -^<p, X(t) = J 1>{s)da, xr(t) = x(|). (3.3) 

We write 

*(x,t) = <p(x)x(t), (3.4) 

* T (x,t) = <p T (x) X T(t), (3.5) 

note that ^ and are positive functions. For generating function F we define 
G(s) = F(l — s) — (1 — s) so in our case 

s i+/3 

G < s > = h? 
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Behavior of the system starting off from a single particle at x is described by 
the function 

(x,r,t) = 1 -Ecxpj-^ (iV;,*(-,r + s))ds| (3.6) 

where N% denotes the empirical measure of the particle system with the initial 
condition Nq = S x . satisfies the equation 

(a;, r,t)= [ %- s (•, r + 1 - s) (1 - v 9 (•, r + 1 - s, s)) - VG (v* (-,r + t - s, s))] (x) ds. 
Jo 

(3-7) 

This equation can be derived using the Feynman-Kac formula in the same way 



as (|15l Lemma 3.4). We also define 

n*(x,r,t)= [ TtsV (-,r + t- s)(x)ds. (3.8) 
Jo 

Since we consider only positive hence (|3.6p and (|3.7p yield 

< v T (x, r, t) < n T (x, r, t), (3.9) 

v T {x,r,t)<l, (3.10) 
In the sequel, for simplicity of notation, we write 

Vt (x, r, t) = Vm T (x, r, t) , (3-11) 

rir {x, r, t) = n-q, T (x, r, t) , (3.12) 

v T (x) = v T {x,0,T) , (3.13) 

n T (x) = n T (x,0,T) (3.14) 

when no confusion arises. 

Fact 3.1. tit (x, T — s, s) — > uniformly in x E K , s E [0, T] as T — > +oo. 

The fact was proved in (|15l . Fact 3.7). From the proof therein we obtain also 
the inequality 

n T < (3.15) 
Following the lines of (flBl Section 3.2.2) we introduce function Vt 

Vt(x,0 = l-Ecxp((A7\ln(l-u T )}) (3.16) 
which satisfies the equation (see (TBI (3.20))) 

y T (x, I) = %v T (x) - V f %_ S G (V T (-, «)) (x) ds. (3.17) 
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A trivial verification using (|3.16[) . (|3.10[) and (|3.17[) provides us with 



< V T (x, I) < %v T (x) , V xm *,l>o- 



(3.18) 



Next we write the Laplace transform of the occupation time fluctuation process 
(|1.5[) for the system N starting from equilibrium distribution 



E e -(* T '*) = e MT)+B{T) 



(3.19) 



where 
A(T) 



(x, T — s)vt (x, T — s,s) + VG (vt {x, T — s, s)) dsdx, 

(3.20) 

r+oo r 

B{T) = V / G(V T (x,t))dxdt. (3.21) 

The derivation of this formula can be found in 

(0, 

Section 3.2.2). 

To show (13.21) we need to calculate limits of A(T) and B(T). Let us notice 
here that exp(A(T)) is the same as the right-hand side of pi (3-7)) so it is the 
Laplace transform of a Poisson-starting system 



A(T) 



V 



1 + 13 

In Section d] we shall prove 

f + oo 



/ p u -r{x)dudsdy 



1+/3 



drdx. 



B(T) 



V 



1+/3 



JR d 



tp(y)pi +s (x)x(s)dsdy) dxdl. (3.22) 



Finally, we are in position to interpret the result. By properties of the Laplace 
transform, the limit of (|3.19jl splits into two independent parts, corresponding to 
A(T) and B(T) respectively. The first was investigated in |7J) and corresponds 
to the process r/ 1 defined by l|2.2j) . 

Let us now investigate the second one. Denote the limit of (|3.22j) by B. It can 
be handled in the following way 

B = / [ ( [ [ Pi+s{x) [ ${y,u)dudsdy) dxdl. 

1 + P Jo jRd \"'M d JO Js J 



We write 



then 



F(x, I) = / / pi+ s {x) / §(y,u)dudsdy, 

JR d JO Js 

B = -r—R / / F{x,l) 1+P dxdl. 

1 + P JO JR d 
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An argument as in (0, Corollary 3.5) implies that the characteristic function of 
is 



exp j-AT f / \F(x,l)\ 1+l3 (l-isgn(F(x,Z))tan|(l + /3)) dxdl\ , 



where A = ^ cos f (1 + /?). 



(3.23) 



1+/3 — 2 < 

Following the lines of reasoning in (7, End of Section 3), we can obtain the 
characteristic function of the finite dimensional distributions (by passing to the 
limit with appropriate sequence approximating $(y, s) — J^j z jVj^t j {s))- It is 
of the form (pQgj) with 



Jo 



Using theorem ( la. Proposition 3.4.2) one can infer easily that rf 1 (recall (|2 .3|) ) 
has the same finite dimensional distributions. 



3.1.2 Tightness 

It has been already mentioned that to prove tightness it suffices to show tightness 
of real-valuedprocesses (Xt,<p) for any (p G 5(K d ). We apply below a scheme 
presented in (0). By |l|, Theorem 12.3) it is enough to show that there exist 
constants \ > and v > such that 

P(| (X T (t 2 ), <p) - {X T {h), <p)\>S)< ^-(h - ti) 1+x (3.24) 

holds for all t 1 ,t 2 G [0,1], ii < t 2 , all T > 1, and all S > 0. A lemma in 
(0, Section 3) shows that each ip G 5(R d ) can be decomposed y> = — <^2, 
yi) V2 G 5(R d ), and (pi,if2 > 0, hence from now on we will assume that <£> > 0. 
Tail probability can be estimated using inequality (0, (3.39)) 

P(| (x T ,ip®<f>) I > 5) < CS 1 (1 - R,e(Ecxp(-i6(x T ,(p®(l^)))d6 (3.25) 

Now an analysis similar to that in (@) shows us (|3.24[) . Indeed, we approximate 
S t2 - 5 tx by G S{R d ) such that x(i) = </.(s)ds fulfils 

0<X<l[ tl ,t 3 ]- 

Notice that /Xt, </? ® 0^ | approximates | (At(£ 2 ), v?) — (JCt^i), </?) |. 
Suppose that we know that 

<S^ 1 (l-Re(m2w(-i0(x T ,tp®<t>})))d9 < ^^-(t 2 - h) 1+x (3.26) 



12 



A passage to the limit using Fatou's lemma implies (I3.24[) . The task of proving 
the last inequality is delegated to Section (|4.3[) . The proof there will be con- 
ducted by estimating characteristic function appearing on the right-hand side 
of (|3.25p . Analogously to (|3.19p we have 

Eexp (-i (x T , ip ® <^>) = exp (A(T) + B(T)) , (3.27) 

where A(T) and B(T) are given by (|3~2T)]) and (pT2"Tj) with v T and V T being the 
complex counterparts of functions v and V from Section 13.1.11 It is easy to 
check that they fulfil equations (|3.16p and 

v T (x,t)= f T t _ s [i<f T C)XT{T-s)(l-VT(;s))-G(v T (;s))}(x)ds 
Jo 

(cf. (32J). 

3.1.3 Auxiliary facts 

Before proceeding to calculations we gather a few additional facts. Let pt denote 
transition density of a-stable motion. We have 

p t (x)=t-«p 1 (xt-i) (3.28) 

And hence 

W 5 = i i(i-«)|| pi ||« )?> ^_. (3.29) 

We will need a few straightforward inequalities 

(a + by+P < 2 1+ P(a 1+ P + b 1+ > 3 ), (3.30) 

(a + b) 1+p - a 1+p - b 1+0 > Pb^a, b > a > 0, (3.31) 
(a + b) 1+p - a 1+p - b 1+p < (1 + /3)aV +/3 " 5 , p < S < 1, a, b > 0, (3.32) 
Moreover, we will use the generalized Minkowski inequality 

II |/|| p < J\\f\\ P ,P>l (3.33) 

and Young's inequality 

H/*5ll,<ll/IUIl5lU, ^ = ^ + ^-1. (3.34) 
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4 Calculations for the proof of Theorem 12.11 



The general schema presented in Section 13.11 left the main technical difficulty of 
the proof, namely convergence of (|3.22p untouched. We decompose B(T) in the 
following way 

B{T) = TTp {B3{T) ~ B2{T) ~ Bl{T)) ' 

where 

r+ao i- 

B 1 (T)= I / (Ttvrix)) 1 ^ -(V T (x,l)) 1+f} dxdl, 

JR d 

+00 







B 2 {T)= / (T inT (x)) 1+ P -(v T (x,l)) 1+ Pdxdl, 

JO JR d 



\l+/3 



dxdl, 



B 3 (T) = / / (Tmrix)) 1 

Jo JR d 

(|3.22[) will be shown once we obtain 

5i(T)^0, (4.1) 

B 2 (T) -» 0, (4.2) 

r+oa r / rl \ 1+ P 

B s {T)^\\ipf+n / / Pl+S (x)x(s)ds) dxdl. (4.3) 

Jo Jm d \Jo J 



The proof of (14. 2|) is omitted since it is similar to the case of (14. 1|) but simpler. 

4.1 Convergence of B 3 

By definition of tit (see (|3.8p ) we obtain 

+00 / T \ 1+13 

B 3 (T)= [ I \Ti I r s ip T (x) X T{s)ds ] dxrfZ. 
Jo Jvi d \ Jo J 

Changing integration variable s — > Ts and using definition of (recall (|3.3|) ) 

B 3 (T) = T°^' f + f ( f T l+T Mx)x(s)ds) dxdl. 
Jo Jvi d \Jo / 

Using the definition of semigroup T (see (|2.1|) ) yields 

B 3 (T)=Ti?- 1 f + f (f f Pl+Ts {x - yMy)x(s)dyds) dxdl. 
Jo Jm d \Jo Jwi d J 

By (pQgj) we have 



B 3 (T) = T^ 1 / + °° / (7 / T-i Pl/T+s (T-i (s - y)) ^(y)x( S )dy<fe 

JO JR d \J0 JR<i v 7 



1+/3 
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and, after obvious substitutions 

B 3 (T)= f + f (f f P i +s (x-T-^y)v(y)x(s)dyd S ) 

JO JR d \Jo JR d K ' J 

This can be written as 



dxdl. 



B 3 (T)= / / (/, * g T ) 1+0 dxdl 
Jo J& d 

where ft(x) = f Q p s+ i(x)x(s)ds and g T (x) = T%(p(xTa). 
Firstly, using the Jensen inequality we easily check 

WfiWltp = II fpi+MWui ^ I [\pi+s(x)) 1+(3 dxds = f \\ Pl+s \\\ + 4ds 



m d Jo 



(by ^M) 



c [\l + S )-^d S < Cl l-iP. 
Jo 

Secondly, using (|3.33p and (|3.29[) we get 



f 1 / f 1 \ 1+0 ( f 1 \ 1+13 

W\\\tt = II I < [J o lbi+.lli+^J = U (i + s)-«^ds) < c 

Combining the last two estimates we get |[/i||Jig < c(l A Z~ a' 3 ). In this way we 
have proved that /; is (1 + /3)-integrable with respect to x and I since — f3 > 1. 

Taking into account the form of qt (informally speaking <?t converges to 
S x ■ \\<p\\) we acquire the £ 1+/3 convergence 

fi*9T-*fr IMI 

which is exactly 



4.2 Convergence of B\ 

We prove (|4.ip . Applying (|3.32p we write 

Bx{T) < B U (T) + B 12 (T), 

where 



/■+OC r 

B n (T)= / {l + p)[TiVT{x)-VT{x,l)}—[Tiv T {x)]—dxdl, 

Jo Js. d 



Bu(T)= / [T lVT {x)-V T {x,l)] 1+0 dxdl. 

Jo JWL d 
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Using (|5T7|l . (j3~Tgj) and (O we get 



/■ + 00 (■ / 

Bi 2 (T) < c / / 



\ 1+/3 

Ti_ s (T s n T (x)) 1+l3 ds ) dxdl 



o 



By definition of (see 

/ r l ( r T 

l (fT{x)du ds dxdl 



£12 (t) < C £ uf 7i_ s ^r s ^ 



Combining with definition of T (|2.1|) and subsituting u — > Tu we can rewrite 

T (i+/3)(#/3-i) / + f (f f pi-.(x-y)(f f P Tu+s(y ~ z)<p(z)dzdu) dyds) dxdl 

JO JR d \JQ JR d \J0 JM d / / 

Application of (|3.28| and substitutions y — > T~iy and x — > Tax yield 

+00 , / rl /■ . . / rl r . . \ H-/J x 1+0 



T c T I I [ I I Pis [T«{x -y)) ( / / Pu+s/t (y ~T « z) tp(z)dzdu) dyds \ dxdl 
Jo JR d \Jo J& d v 7 Vo JR d v 7 / 



where C T = -(1 + (3) 2 + & (1 + /3 2 + /3 3 ). Applying and substituting 

s — > Ts and / — > TZ we obtain 

7 J / [iff pi- s (x-y)fj J p u+s (y -T~^zj (p(z)dzdiij dyds] dxdl 



10 Js. d Wo 
Let us denote 

• z ) <p{z)dzdu 



h{s,y) = ( / / Pu+s (y — T a 

\Jo JR d K 



1 x 1+/3 

*T * / Pu+sdu I (y), 



11 



where $t(x) = T<*ip(T"x). We obtain 

B 12 (T) < C T / ff(0<fl, 
Jo 

where 



(4.4) 



#(0 = 11 / Pl-s*Ms,-)ds\\^. 
Applying (|3.33p we get 



(i 
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Utilizing Youngs 's inequality (|3.34|) we write 

H(l)< ^||p;- s || 1+/ 3||/t(s,-)||i^ ■ 

By (|3.29|) and the definition of fx we obtain 

/ i 1 \ 1+13 

H(l)< lj(l-s)-^U\\$ T *^ p u+s du\\ll^ds) . 

Using (|3.34p once again we have 



h(i)< iy o a-*r^n*Tiir p ii y o pu + ,^ii^j 

Hence by the (1333]) and (f3T29|) 

H(l)<c(j (I - s)~ (J (u + s)~^+^du) ds 



1+13 



For intermediate dimensions — -tt^tf < 1, hence the inner integral can be esti- 

a 1+p ' & 

mated by a constant independent of s so 

/ i \ 

#(0 <clj(l- sy^mdsj 

Assume Z < 1. Then 

< c. (4.5) 
Now we derive estimation that works for "large" Vs. By (|3.33[) we have 

H(l)<U\\ Pl _ s *f T (;s)\\ 1+0 d S ) . 
Young's inequality (|3.34|) yields 

H{i)< ^'||pi-IIi + /9II/t(-,«)|| 1 ^ • 

We estimate C 1 norm of f t using (|3.34[) 

n/ T (., S )ii 1 = ii$ r * f 1 Pu+S du\\ix ^w^wi^w f Pu+s du\\\ + + % 

Jo JO 
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We use the trivial fact that ||$t||i = c an d (|3.34[) (p, q G [1,1 + /?] are yet to be 
specified) 

< ll„ lli+z 3 !! /~<n W^lli+z 3 



< C||p s * / VuH\\X% < \\Vst+ P \\ \ Pudu\\\ . 

Jo Jo 

We have 

II jf Puduf+^K (J \\ Pu \\ q du\ < (J u^-^du) <c, 
since it is obvious that for any q G [1, 1 + 0\ we have ^(^ — 1) > — 1- And finally 

HM-.-)lli<«* ( *- l)(l+/,) - 



One can adjust p to make exponent A = — f^ — lj(l+/3) arbitrary near —1 



d II 

a \p 

(because if p = 1 + /3 then /I < —1). Hence going back to estimation of H(l) 
H{1) <c( [ (I- s) B s A ds] 



where B = — ( y^-g — 1 ) • Substituting s — > we obt 



am 



1+/3 , „i x 1+/3 



Notice that B > — 1 and 

(A + S + 1)(1 + /3) = (A+ 1)(1 +0) + (--(3 



~a@ ^ — ^ an d ^ + 1 can be made arbitrarily near 0, hence 

H(l) < l w (4.6) 



where W < — 1. Combining estimates (|4.5p and (|4.6[) for if we can conclude 
that the integral in (14. 4| is finite and £?i2 — » 0. 



4.3 Tightness calculations 

Following the scheme in Section [3.1.2l it remains to prove (J3T26J) . Slightly abusing 
notation we will use an additional argument 9 to indicate that a function is 
computed for instead of $ (eg. A(T,6)). In this section we deal with 
complex functions defined in Section f3.1.21 which are not to be confused with 
functions in sections devoted to space-time convergence. 
From equation (|3.26p we have to estimate 



L = 1 — Re [ Ecxp -i6 ( X T , tp ® <j) 
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We use Q3.27I) and then estimate A(T, 6) in the same way as in the proof of 
tightness in (j7|), thus 



L < E 
where 



3> £ 

1 - exp (-i6 (x T , ip ® </>\) < |i + u + iu| 



I = i0 / (/3t(2;)xt(7 1 — s)i>T,e(a;, s)dxds 

JR d JO 

II '= T-^—J / / V^{x,S)dxds 
1 + P JV JO 

III = V / (Vr,e(x,i)) 1+/3 dxrft 

JO JR d 

The terms / and II are the same as in ([j]). Hence we have only to deal with 
\III\- Before that we show an estimation (which holds for T large enough). 
Firstly, recall the definition of Vt (|3. 16[) 

\V T ,o{x,t)\ = |1 -Ee^'^^^l <E|l-e (Ar * X ' lnu,r . 8> | 

where t 

WT,e(x,r,t) = Eexp< -i6 / (JV*, * (■, r + s)) ds 



We know that |wr,e| < 1 which implies jlnwT.el < and consequently e^^' 1 "™^^ < 
1 . Finally, if \z \ < 1 we can use inequality 1 1 — e z \ < 2 \ z \ . Hence 

\V T ,e( x ^)\ <2E|(JVf,lnt«r, fl >| < 2E{N? , \ \nw T ,e\) = 2T t \ \nw T \ <2T t n T . e . 

(note that riT,e is a real function). 
Therefore we have to estimate 



p + ca p 

\Bt,b\ < / / {Tin T ,e{x)) 1+l3 dxdl 
Jo Ju d 



Let us notice that the integral is the same as B3 from in Section [4j For T large 
enough we have 



1+/3 



\B T ,e\ < C\\eip\\ 1+/3 / / ( p l+s (x)x(s)ds) dxdl 
Jo JR d \Jo / 

According to the argument in Section [3. 1.21 we choose \ — l[ti,t 2 ] 

,t 2 x 



|S T ,e| < 
Denote A := <2 — £i 



/ / ( / Pi+s ( x ) ds ) dxdl 
Jo Jni d VJti / 



Bt,8 < / / I / Pl+s ( x ) ds rfxcW 

JO JR d \J0 / 
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After obvious substitutions and using (|3.28[) we obtain 

\BtM < A< 2 +«A^^ / / Pl+S (x)ds) dxdl. 

JO JR d \Jo / 

It is easy to check that 2 + (3 — ~(3 > 1 reasoning along the lines of the proof in 
(0) completes the proof. 
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